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ABSTRACT 

The importance of the wall effect on packed beds in the adsorption and desorption of carbon dioxide, 
nitrogen, and water on molecular sieve 5 A of 0.127 cm in radius is examined experimentally and with one- 
dimensional computer simulations. Experimental results are presented for a 22.5-cm long by 4.5-cm 
diameter cylindrical column with concentration measurements taken at various radial locations. The set of 
partial differential equations are solved using finite differences and Newmans’ s method. Comparison of test 
data with the axial-dispersed, non-isothermal, linear driving force model suggests that a two-dimensional 
model (submitted to Separatuin Science and Tchnology) is required for accurate simulation of the average 
column breakthrough concentration. Additional comparisons of test data with the model provided 
information on the interactive effects of carrier gas coadsorption with C0 2 , as well as C0 2 -H 2 0 interactions. 
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Introduction 


In a packed bed, the presence of the confining walls can have a significant influence on the axial flow 
and therefore heat and mass transfer in the packed bed. The viscous shear force along the wall and the 
increase in porosity near the wall can have opposing effects on overall permeability. For values of tube to 
particle diameter ratio, R y less than about 7, the permeability has been shown to increase above that of an 
infinite porous medium (1,2,3). Between R values of 7 and 30, variations in bed packing may determine the 
governing influence. For R values above 30, the influence is generally considered insignificant (3), although 
it was found to be important for an R value greater than 100 in an adsorption study by Tobis and Vortmeyer 
(4). 

Wall effects on packed bed processes can be significant in heat and mass transfer driven processes, since 
the value of R is frequently within this sensitivity range. Examples are thermal swing adsorption beds where 
a heater matrix results in small rectangular columns (CDRA refs), wall-cooled catalytic reactors (5) and 
packed bed energy storage units (6). 

In this work, we examine the influence of wall effects on mass and heat transfer during adsorption and 
desorption in a packed column. The experimental results for carbon dioxide and water adsorption in a 
nitrogen carrier gas provides information on the sensitivity of the adsorption to wall effects in a column 
with an intermediate tube to particle diameter ratio, R y of 15. Adsorption breakthrough data is presented for 
various radial locations at the bed exit, which can be used as a measure of permeability or channeling along 
the column wall. 

A computer model was developed based on the finite-differencing numerical technique. The model 
simulates gas adsorption and desorption in a flow-through bed. Single and multicomponent adsorption; 
thermal balances for the pellet, gas and canister, and the momentum equation are included. A two- 
dimensional (radial and axial) model was also developed to predict the effect of flow channeling on 
adsorption and desorption and is described in a later paper (7). 

Experimental Apparatus and Results Adsorption Test System 

The test flow rate was scaled to give gas velocities similar to that in the proposed flight beds for the 
International Space Station Carbon Dioxide Removal Assembly. Instrumentation of the packed column 
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includes temperature probes and sampling tubes for measurements at sorbent material endpoints and one 
intermediate point. To attain approximately adiabatic bed conditions, the column was first covered with 
Mansfield Q-fiber felt insulation, then wrapped with a thermal blanket of Mansfield Min-K material. For 
adsorption runs, an additional jacket fabricated of STS External Tank insulating foam was attached. The 
insulation, approximately 3 inches thick in all, was used with satisfactory results. 

Instrumentation was provided for continuous measurement of packed column outlet C0 2 and H 2 0 
concentration. The location of sensors and other equipment comprising the adsorbing apparatus are shown 
in Figure 1. The column bypass is used prior to the test start to ensure stable column inlet conditions. 

The gas chromatograph, a Shimadzu GC-14A with CR601 integrator, was used to sample at three 
sample port locations to determine gas constituent volumetric fractions during the adsorption runs. The 
probe depth of the sampling tube at the column exit was adjusted to obtain a radial profile of the exit 
concentration during a series of identical C02 adsorption tests. 

Column Dynamics Test Bed 

A small packed column with approximate dimensions of 2 inches in diameter and 20 inches long (Figure 
2) located at Marshall Space Flight Center was used. The column may be packed with up to 20 inches of 
sorbent. A 10-inch packing was used in this study for quicker results and reduced thermal end effects. 
Placing 4.75 inches of glass beads at the two ends of the column eliminated the end effects. 

Procedure 

Experiments were performed on the insulated fixed-bed rig at MSFC. The sorbent is a 5A zeolite. 
Nitrogen is the carrier gas. Experiments begin with the column at ambient room temperature. The column 
bypass shown in Figure 1 was used to obtain intended inlet conditions before exposure of the column to the 
inlet gas. The GC sampling location was switched during the test to follow the sorbate mass transfer wave 
as it proceeded down the bed. 


Results and Discussion 

Figure 3 shows the water breakthrough at the column midpoint and outlet. Also shown is the H z O partial 
pressure of gas mixed by passing through glass beads downstream of the sorbent material. Note the mixed 
gas breaks through before the gas at the centerline, indicating that channeling has a significant effect on the 
process efficiency for the two-inch diameter column. 
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To provide a measure of the radial variation in axial flow, adsorption tests of C0 2 and N 2 were 
conducted repeatedly with the exit sampling tube in different radial positions. The test results, shown in 
Figure 4, show a clear relationship between concentration and radial position. As observed for H 2 0 
breakthrough results, these results indicate greater permeability near the column wall. Since the heater core 
of the 4BMS (Four Bed Molecular Sievs) sorbent bed consists of channels roughly one-half inch in 
diameter, the channeling effect on C0 2 removal will likely be significant. The results of the two- 
dimensional model presented in part II also confirm that channeling has a significant effect on this 
adsorption process. 

Mathematical Model 
Gas/Solid Equilibrium 

The model discussed here was developed using equilibrium data provided by the W.R Grace & 
Company (8). Nitrogen adsorption was represented with the Langmuir expression: 
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Carbon dioxide and water adsorptions were correlated by the Langmuir-Freundlich equation (9): 
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Two approaches were used to solve for the amounts of gases adsorbed in the solid phase in equilibrium 
with the gas phase concentrations for multiple components. For simulations with N 2 and C0 2 , the solution 
was obtained by an iterative procedure with the ideal adsorption solution theory (IAST) of Myers and 
Prasusnitz, (10). For the case of the Langmuir-Freundlich single component isotherm, one gets: 
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Knowing the total pressure and the solid temperature, a value for the spreading pressure k is 
estimated and solution will be obtained by iteration until Eq. 3 is satisfied. For the simulations with N 2 , 
H 2 0, and C0 2 the approach based on the Langmuir-Freundlich equation for multicomponent systems was 
found to give a more accurate results than IAST model based on the non-ideality of co 2 /h 2 o, a paper is in the 
process for submision to discuss in details the use of IAST and Langmuir-Freundlich on the adsorption of 
C0 2 /N 2 and H 2 0/C0 2 /N 2 on the 5a : 
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The favorable agreement of simulation and test results using these correlations indicates low sensitivity to 
limitations of these correlations for components with differing adsorptivities. However, work continues 
towards a comprehensive gas/solid equilibrium model. For example, at low humidity levels, IAST 
accurately predicts multicomponent adsorption equilibria for C0 2 and H 2 0. At higher humidity levels, C0 2 
adsorption is under predicted, as discussed by Finn (11). 


Mass Balance Equation 

In the bulk stream of the gas within the bed, the material balance for the adsorbate concentration is 


3 Ci d^Cj duCj 1 -edqi 

3 1 l 9^2 dx £ dt 


Boundary and initial conditions: 

at t < 0, Ci = Ci q for 0<x<L 
at t < 0, qi -qi y Q for 0<x<L 

att>0 , Ci = Ci i for x~0 ^ 

dC* 

att >0,-^- = 0 for x= L 
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The axial diffusion, D h was calculated using an equation from Edward and Richardson equation (12). 

Assuming the ideal gas law C i = P / RT and knowing ^ V. = 1 , the above equation can be recast 
into an overall mass balance equation, 
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This equation was used to compute axial velocity in the bed. 

Gas Phase Energy Equation 

The change of gas temperature with respect to time is due to heat flux from the solid to the gas plus 
convection of heat due to the fluid flow, as shown by the following equation: 
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Boundary and initial conditions: 


at t <0,Tg = Tg q for 0<x<L 
at t > 0, Tg = T( for x-0 
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Note that T w is calculated in Eq 12. 


Solid Phase Energy Equation 

The energy equation for the solid phase includes the term for heat flux from the solid phase to the gas 
phase plus heat generation due to adsorption. 
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Column Wall Energy Equation 


The wall temperature T w is given by: 
clT 

P w C pw "af" = 2nR < h ^ T x ~T w )-2nR 0 h 0 (T w -T 0 ) (12) 

Initial condition: 

att<0,T w =T Wt0 (13) 

Axial conduction is neglected since the area of heat transfer from the fluid to the wall is an order of 
magnitude larger than the area in the axial direction. This is analogous to heat conduction in a slab. 


Momentum Equation 

The Ergun equation (13,14) is used to estimate the pressure drop, 

dP /i ~ 2 

= u — pCu 

dxK 


(14) 


where C is the inertial coefficient. The empirical coefficients K and C are given by relations developed 
by Ergun for flow in a packed bed of spherical particle; 
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Using gas velocity as calculated in Eq. 7, the total pressure was found with Eq. 14. 

Solid Phase Transport Equation 

Mass transfer of solute from the bulk gas to sorbed state is driven by the difference in the actual 
adsorbed quantity versus the quantity that would be adsorbed at equilibrium conditions. In general, the mass 
transfer mechanism of an adsorption process includes four steps: fluid-film diffusion, pore diffusion, 
adsorption rate, and surface diffusion. Adsorption rate can be neglected since it is much greater than the 
diffusion rates as discussed by Yang, (15). 

Zeolite sorbents consist of crystals, in the size range of 1-9 pm, which are pelletized with a small 
amount of binder. Macropore (spaces between the crystals) and micropore (intracrystalline) diffusion must 
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in general be considered. However, for molecules in the size range of C0 2 adsorbing onto 5A the diffusion 
rate inside the pellet has been shown to be controlled by intracrystalline diffusion (16). 

Assuming the intracrystalline diffusion governs the overall mass transfer, the linear driving force 
(LDF) model based on Glueckauf, (17) may be used: 


A 


= k ef o s (q i -qi) 


(17) 


where k ef is obtained by experimental procedure and a s is the interfacial surface area. The justification of 
assuming a linear diving force to model the adsorbed concentration in the solid phase has been well 
established by other researchers such as, (16,21 ,22,23) to name a few. 

Gas to pellet heat transfer coefficients obtained from the Petrovic and Thodos (18), for C0 2 , H 2 0, 
and N 2 x were .0875, .1002, and .0844 cm/sec, respectively. The empirically obtained mass transfer 
coefficients using the LDF model for C0 2 , H 2 0, and N 2 were 1.27X10 -4 , 2.625xl0' 5 , and 3.75xl0‘ 4 cm/sec, 
respectively. Since the overall mass transfer rate is three orders of magnitudes smaller than the calculated 
mass transfer rate from the gas phase to the pellet surface, the latter was neglected. 

Numerical Solution 

The solution of this system must be obtained numerically. For an n component mixture, the numerical 
model will require the solution of several coupled differential equations: n-1 mass balance equations, n rate 
equations and equilibrium isotherms, the overall mass balance equation, the momentum equation, and the 
heat balance equation each for the fluid phase, gas phase, and column wall. 

In this work, the PDE’s were discretized using the finite difference method. The first order 
approximation was used for the time dimension. First order boundary nodes and second order internal nodes 
were used for spacial coordinates. The set of discretized finite difference equations was solved 
simultaneously by an implicit method. The method of Newman (19,20) was found to be effective for the 
steep gradient of mass concentration along the bed length, Newman, 1967and 1968. 

50 nodes represented the column length. The time step was increased from 3.6 to 60 seconds as the 

i • _ /i+l n n+1 . 

solution progressed. Convergence was assumed when (C +C )/C was less than 1.0x1 0' 4 for each grid 
point. Here, C is the gas phase concentration, n indicates the previous time step, and n + 1 the current time 
step. 
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Results 


In this section we present modeling results which illustrate the importance of accurately simulating 
critical adsorption processes, and comparisons of model results with test data. Simulations were performed 
on an ALPHA VAX computer. Experiments were performed on the insulated fixed-bed rig at MSFC. The 
sorbent is a 5A zeolite. Nitrogen is the carrier gas. Experiments begin with the column at ambient room 
temperature. 


Thermal Effects Modeling 

Heat transfer coefficients are found by comparison of the thermal model with heating of the sub-scale 
column with an inert gas, shown in Figure 5. 

Nitrogen was heated to a temperature of 350°F and passed through an initially cold column. The 
markers are the measured temperatures at the inlet, midpoint, and exit of the column. The lines are the 
results from simulation. Resultant heat transfer coefficients are 2.5 BTU/hr-ft 2 from the fluid stream to the 
canister wall, and 0.25 BTU/hr-ft 2 from the canister wall to the atmosphere. These values are used in the 
simulations described below. 


Nitrogen Coadsorption Effects Modeling 

The importance of not neglecting the nitrogen co-adsorption with C0 2 is shown in Figure 6. C0 2 in 
nitrogen at 6.2 mmHg was passed through an initially clean zeolite column. The gas at the column 
centerline was analyzed periodically at the column inlet, midpoint, and outlet. These results are compared 
with simulations that either included (solid lines) or neglected (dashed lines) the effect of nitrogen. Nitrogen 
clearly has a noticeable effect; however the simulations that included N 2 adsorption over-predicted the 
effect slightly. 

The thermal sensitivity to N 2 co-adsorption is shown in Figure 7. The agreement is much better when the 
heat of adsorption of nitrogen is included. A decrease in temperature occurs at 0.15 hours in both the test 
data and the simulation, due to N 2 displacement by C0 2 and N 2 reduces the bed capacity and lower overall 
bed temperatures through desorption. 
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Water Coadsorption Effects Modeling 

Verification of the model for C0 2 /H 2 0/N 2 co-adsorption is shown in figures 8 to 10. These three figures 
illustrate testing and simulation runs for adsorption of water at 6.3 mmHg and C0 2 at 2.89 mmHg in a 
carrier gas of nitrogen. All figures show test data as markers and simulation data as lines. 

Figure 8 illustrates the roll-up phenomenon as adsorbed C0 2 is driven off by water, and the C0 2 partial 
pressure rises above the inlet level temporarily. The effect should be larger at the column outlet than 
midpoint as shown by test data. 

Figure 9 shows the water breakthrough at the column midpoint and outlet. Also shown is the H 2 0 partial 
pressure of gas mixed by passing through glass beads downstream of the sorbent material. Note the mixed 
gas breaks through before the gas at the centerline, indicating that channeling is significant along the walls 
of the two-inch diameter column. Since the heater core of the 4BMS sorbent bed consists of channels 
roughly one-half inch in diameter, the channeling effect on CO2 removal will be significant. The results of 
the two-dimensional model (not shown here) confirm that channeling has a significant effect on this 
adsorption process. Finally, Figure 10 shows the temperatures of the bed at midpoint and outlet. 

Bed Regeneration 

Desorption of C0 2 test results in comparison with model are shown by the solid lines in Figures 1 1 and 
12. After the bed was saturated with CO2, the regeneration process was started by using N2 as the purge 
gas. As it shown the effluent concentration of CO2 has sharp drop at first few minutes and the slope of the 
breakthrough flats out as time goes by. The initial drop in temperature is due to heat of desorption and 
finally reaches the inlet condition when there is no depletion of CO2 from the bed. The same mass transfer 
coefficient of .017 was used. The model predicts both temperature and breakthrough fairly well. The I AST 
was used to predict the mixtureisotherm of CO2/N2 

Desorption of H 2 0 test results in comparison with the model are shown by the solid lines in Figures 13 
and 14. The results of the partial pressure of H2O and temperature profile of this 1-D desorption model are 
not in good agreement with the test data. The model shows a fast reduction of H2O partial pressure in gas 
phase. A mass transfer coefficient as large as .04 ft/hr, in contrast with .0035 ft/hr in the case of H2O 
adsorption, was used. Generally, there shouldn't be such a large difference between the two coefficients. 
Even with this large mass transfer coefficient the desorption of H2O from the bed was insufficient to 
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increase the H2O partial pressure in the gas phase. In contrast with adsorption process, any small 
discrepancy of H2O partial pressure with test data will remain as a error throughout the completion of the 
test. In adsorption any small error at any point in the bed, if it is caused by the isotherm at some partial 
pressure of the feed will be eliminated at a later time because of the correct value of isotherm at a larger 
partial pressure of the feed. This can be seen from the early breakthrough observed by other researcher. The 
obvious reason is that the equilibrium isotherm at low partial pressures are being underestimated. It is also 
possible that the equilibrium isotherm of H2O on 5 A material shows hysteresis so that the apparent 
equilibrium pressure observed in adsorption and desorption experiment are different. The concentration of a 
key component, CO2, is affected by the presence of the non-key component, N2, in CO2/N2 adsorption. 
CO2 effluent concentration overshoots its inlet concentration because of H2O displacement (more easily 
adsorbed) component. The height of this roll-up is increased with the inlet concentration of H2O 
component. The most significant contribution to the difference in model and experimental results of H 2 0 
desorption pretend to that fact that the duration of an H 2 0 adsorption run takes about fifteen hours for 
completion of a test. Also, the desorption run duration takes ten to twelve hours. During this long duration 
the temperature of location where the test took place varied about 10-15 degree °F, during the night. This in 
turn affects the saturated air that was used to saturate the column. Therefore it is accurate to conclude the 
amount of H 2 0 adsorbed on the bed is less than what is assumed in model. 

On the basis of the data presented here, and other comparisons between desorption test data and 
simulation results not yet published, this computer model meets its primary objective - achieving predictive 
capability. Enhancements to the model as discussed should increase its accuracy. Efforts are continuing to 
develop on the integrated 4 BMS simulation, equilibrium isotherms, heat and mass transfer coefficients, and 
verification data. 

Conclusions 

Based on analytical and experimental investigation of convective flows in porous media the 
following conclusions are drawn: 

• The experimental results from the laboratory scale-fixed bed adsorber are quantitatively consistent with 
the one-dimensional model at the column center. The average concentration of cross sectional bed 
obtained by test result deviates from the column center concentration appreciably. This indicates the 
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strong effects of porosity variation along the radial direction of column bed on the temperature, 
concentration, and velocity field. Results from the model were encouraging and contributed to the 
decision to model the dynamic behavior of the column in two-dimensions. 

• A linear driving force mass transfer model provides a reasonable fit to experimental adsorption data. 

• The concentration of a key component, CO2, is affected by the presence of the non-key component, N2, 
in CO2/N2 adsorption. CO2 effluent concentration overshoots its inlet concentration because of H2O 
displacement (more easily adsorbed) component. The height of this roll-up is increased with the inlet 
concentration of H2O component. 


NOMENCLATURE 

A Surface area of pellets per unit volume of pellet ft 2 /ft^ 

A Surface area ft 2 

B Langmuir constant 

C Constant in Darcy equation 

C Gas stream concentration, lb mole/ft^ 

c ip Gas phase concentration of i^ 1 component in the pores of particles, lb mole/ft^ 

Cip Volume average pore Gas phase concentration of i^ 1 component i, lb mole/ft^ 

C- Gas phase concentration of i^ 1 component at the surface of the solid phase lb mole/ft^ 

Ci q Gas phase concentration of 'fi 1 component at boundary or initial lb mole/ft^ 

Cpg Heat capacity of gas phase, Btu/lbm-R 
Cps Heat capacity of solid particle, Btu/lbm-R 
Cpw Heat capacity of column wall , Btu/lbm-R 

D diffiisivity, ft 2 /hr 

D/ Axial diffusion, ft 2 /hr 

H 0 Effective Heat transfer coefficient for column insulation, Btu/ft 2 -hr 
H w Heat transfer coefficient between the gas stream and the column wall, Btu/ft 2 -hr 
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Hs Heat transfer coefficient between the gas stream and the sorbent , Btu/ft^-hr 

A H Heat of Adsorption , Btu/lb mole 

K Constant in Darcy equation 

Kf axial conductivity of fluid flow, Btu/ft-hr-R 

K s Solid thermal conductivity, Btu/ft-hr-R 

A/ ; Molecular weight of adsorbate i, lb/lb mole 

N Number of component 

P Total pressure, mmHg or lbf/ft 2 

Pi Partial pressure of component i, mmHg or lbf/ft 2 

Q Amount adsorbed in the solid, lb moles/ft^ of solid 
* 

Solid phase concentration of i^ 1 component in equilibrium with gas phase, lbs moles/ft^ of solid 
q m i Langmuir constant 

R Ideal gas constant 555 mm Hg ft 3 /lb mol R 

Ri Inside wall diameter of column, ft 

R 0 Outside wall diameter of column, ft 

Rp Particle radius, ft 

T Time , hr 

T Temperature R 

T 0 Ambient temperature, R 

Tg gas temperature, R 

T w Wall temperature, R 

T s Solid temperature, R 

U Interstitial velocity, ft/hr 

Greek Letters 

£ External bed void volume 

Ppp Density of gas phase, lb mole/ft^ 
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p s Density of solid phase t lbs/ft^ 

Pw Density of column wall, lb mole/ft^ 
K Constant or spreading pressure 

AH Heat of adsorption, BTU/lb of solid 

Subscripts 


i 

i 1 * 1 component 

e 

Effective 

0 

Outside, initial 

P8 

Gas phase 

ps 

Solid phase 

s 

Surface 

t 

total 

w 

Wall 


Superscripts 

Average value 

* Equilibrium value 
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Figure 2: Column Sensor and Sampling Tube Location 
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Figure 5: Thermal Model Validation 
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Figure 6: Effect of Nitrogren Coadsorption on C02 breakthrough 
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Figure 8: C02 breakthrough for C02/H20/N2 Coadsorption 
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Figure 9: H20 breakthrough for C02/H20/N2 Coadsorption 


25 








Figure 1 1 : C02 Depletion for C02/N2 desorption 
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Figure 13: H20 Depletion for H20/N2 desorption 
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Figure 14: Heat Transfer Effects of H20depletion for H20/N2 desorption 
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